Noname manuscript No. 

(will be inserted by the editor) 



Analysis of dislocations in quasicrystals composed of 
self-assembled nanoparticles 

Liron Korkidi • Kobi Barkan • Ron Lifshitz 



the date of receipt and acceptance should be inserted later 



Abstract We analyze transmission electron microscopy (TEM) images of self- 
assembled quasicrystals, composed of binary systems of nanoparticles. We use an 
automated procedure that identifies the positions of dislocations and determines 
their topological character. To achieve this we decompose the quasicrystal into its 
individual density modes, or Fourier components, and identify their topological 
winding numbers for every dislocation. This procedure associates a Burgers func- 
tion with each dislocation, from which we extract the components of the Burgers 
vector after choosing a basis. The Burgers vectors that we see in the experimental 
images are all of lowest order, containing only O's and I's as their components. We 
argue that the density of the different types of Burgers vectors depends on their 
energetic cost. 



1 Dislocations in self-assembled soft-matter quasicrystals 



Self-assembled soft-matter quasicrystals have been observed in recent years in a 
wide variety of different systems, in all cases but one with dodecagonal (12-fold) 
point-group symmetry. First discovered by Zeng et al. (2004) in liquid crystals 



made of micelle- forming dendrimers, self- assembled soft-matter quasicrystals have 
since appeared in other systems such as ABC-star polymers (Hayashida et al. 
2007), in binary systems of nanoparticles ( Talapin et al.|20Q9 ), in block co-polymer 
micelles (Fischer et al. 2011), and in mesoporous sihca ( |Xiao et al.|[2QT2l ). These 
newly-realized systems not only provide exciting platforms for the fundamental 
study of the physics of quasicrystals (Barkan et al. 2011), they also hold the 



promise for new and exciting applications, especially in the field of photonics. An 
overview of soft matter quasicrystals, including many references relevant to these 
systems, is given by Lifshitz and Diamant (2007) as well as by Ungar et al. ( 2005[ 
20TT| and |Dotera| ( |2011| |2012f 

Here we concentrate on the systems of nanoparticles studied by [Talapin et al.| 



(2009), consisting typically of two types of particles, such as PbS, Au, Fe203, 
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and Pd, with different diameters. These binary systems of particles, when placed 
in solution, self-assemble into structures with long-range order, including 12-fold 
symmetric quasicrystals. The dimensions of the particles — typically a few nanome- 
ters in diameter — are such that they can be imaged directly using a transmission 
electron microscope (TEM). This allows one to study effects that are inaccessi- 
ble with atomic-scale quasicrystals. Here we present a quantitative analysis of the 
dislocations that are naturally formed in these quasicrystals as they self-assemble. 

In periodic crystals in d-dimensions one can usually identify the position of a 
dislocation rather easily by the termination of a plane of atoms in three dimensions, 
or a line of atoms in two dimensions. One then chooses a basis for the periodic 
lattice; encircles each dislocation with a Burgers loop, or a Burgers circuit, of basis 
vectors; and counts the accumulated difference between the number of steps taken 
forward and backward in the direction of each of the d basis vectors. The d integers 
thus obtained define the Burgers vector which encodes the topological character 
of the dislocation. A similar real-space procedure can be used on a quasiperiodic 
crystal by overlaying it with a quasiperiodic tiling of rank D > d (for a definition, 
see below), yielding a D-component Burgers vector. A tiling-based analysis of 
binary systems of nanoparticles was indeed recently carried out by |Bodnarchuk| 



et al. (2011). Here we propose an alternative approach for analyzing dislocations 



in Fourier space that we believe is useful when dealing with aperiodic crystals. 



2 Density modes, winding numbers, and the Burgers function 

Let us describe the density of nanoparticles in a self-assembled crystal by a function 
p(r). The Fourier expansion of such a function is given by 

p(r) = ^p(k)e^'^-, (1) 

where the (reciprocal) lattice L is a finitely generated Z-module, which means that 
it can be expressed as the set of all integral linear combinations of a finite number D 
of d-dimensional wave vectors, b^^^ , • • • , b'^^^^ . In the special case where the smallest 
possible D, called the rank of the crystal, is equal to the physical dimension d, the 
crystal is periodic. More generally, for quasiperiodic crystals D > d, and we refer to 
all quasiperiodic crystals that are not periodic as quasicrystals (see Lifshitz||2QQ3| 



2007). 



As explained elsewhere (Lifshitz||2011 ), the topological nature of a dislocation 



is related to the fact that it cannot be made to disappear by local structural 
changes. For this to be the case, as one follows a loop around the position of a 
dislocation and returns to the point of origin, one sees a crystal that is everywhere 
only-slightly distorted from the perfectly ordered state, except near the core of 
the dislocation. In particular, the complex amplitudes p(k) of the density modes 
maintain their magnitudes along the loop, each accumulating at most a phase, 
which upon return to the point of origin must be an integer multiple of 27r. The 
collection of all such integers, or so-called winding numbers, for a given dislocation 
defines a linear function MbO^) from the lattice L to the set of integers Z, which 
we call the Burgers function. 

The Burgers function of a given dislocation associates a particular winding 
number A/j5(k) with every wave vector k G L. Because this function is linear, after 
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Fig. 1: (a) A section of a TEM image of a dodecagonal quasicrystal, self-assembled 
from 11.2nm PbS and 5.2nm Au showing several dislocations (courtesy of Dmitri 
Talapin 2012). (b) Fourier transform of the TEM image in (a), with the central 



peak blocked. A pair of Bragg peaks, associated with one of the basis vectors and 
its negative, is marked in red. (c) The corresponding section of the inverse Fourier 
transform of the Bragg peaks marked in (b) with red circles marking the positions 
of four dislocations. 



choosing a basis {b^*^} for the lattice, it is uniquely specified by a set of only D 
integers = A/j[3(b^*^), forming the Burgers vector (ni, . . . ,n]j). Thus, 



D 



D 



D 



Vk: 



^a,b^^^ G L : A/-^(k) = ^ a,A/-^(b^^)) = ^^a.n, 



(2) 



i=l 



where G Z. This implies that in order to fully characterize a dislocation in an 
experimental image it suffices to isolate the D density modes associated with a 
chosen basis, and obtain their corresponding winding numbers. This is the basis 
of the approach presented below for analyzing dislocations (for more detail, see 
Barak and Lifshitz||2006[ [Freedman et al.||20Q6[ |2Q07[ ). 



3 Analysis of the dislocations in a quasicrystal of nanoparticles 

We begin with a high-resolution TEM image of one of the dodecagonal quasicrys- 
tals grown by Talapin ( 201 2 ( ), a section of which is shown in Fig.[TJa). This partic- 
ular quasicrystal is self- assembled from 11.2nm PbS and 5.2nm Au nanoparticles, 
and contains a distribution of dislocations that are formed naturally during the 
self assembly. We Fourier transform the TEM image, to obtain the diffraction im- 
age shown in Fig. [ijb), and then choose four of the Bragg peaks in the 12-fold 
ring containing the strongest reflections as a basis h'^^\ . . . , b^'^'^ for the reciprocal 
lattice. These are labeled in the schematic representation of the lattice in Fig. [2] 
For each of the four pairs of Bragg peaks, associated with the chosen basis 
vectors and their negatives, we then carry out the following procedure: 

1. We filter out small regions around the two opposite Bragg peaks, as indicated 
by a pair of red circles in Fig. [ijb) for the case of the density mode associated 
with the wave vectors ±b^^^ 
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Fig. 2: Schematic representation of the three strongest rings in the Fourier trans- 
form of our dodecagonal quasicrystal. The inner ring is the strongest, containing 
the four basis vectors b^^^ . . . h^^\ The second strongest ring is the outer one, ob- 
tained from ah the sums of two adjacent vectors in the inner ring, as indicated 
by sohd arrows. The third strongest ring hes in between, obtained from all sums 
of vector pairs in the inner ring that are separated by 90 degrees, as indicated by 
dashed arrows. 



2. We inverse Fourier transform the filtered regions resulting in a real-space image 
of a single density mode. Dislocations appear as discontinuities in the stripes. 
We use a routine that identifies all the discontinuities and marks their positions, 
as shown in Fig. [ijc) for this density mode. 

3. For each dislocation, a second routine then extracts the i^^ component rii = 
A/'j[3(b^*^) of the Burgers vector. This is done by enclosing a counter-clockwise 
loop around its position and calculating the accumulated phase in units of 27r. 
Practically what we do is count the number of stripes crossed moving in the 
direction of the wave vector b^*'^ on one side of the dislocation, and subtract the 
number of stripes crossed moving against the direction of b^*^ while returning 
on the other side, as demonstrated in Fig. [3] 

Finally, we verify the correctness of the calculation by extracting the values A/j[3(k) 
for additional wave vectors k and checking that they satisfy the linearity require- 
ment given by Eq. 



4 Results and discussion 



We typically find a density of a few dozen dislocations per /xm^ in the nanoparticle 
quasicrystals of Talapin (2012). All of these dislocations are of lowest order in 
the sense that rii = 0, 1 or —1. To understand the topological nature of these 
dislocations it is useful to classify them by dividing the four basis vectors into 
two hexagonal subsets — {b^^\b^^^} and {b^^\b^^^} (see Fig. [i]). By doing so we 
find that the density of dislocations with non-zero components in only one of the 
subsets, which we call single-subset dislocations, is five times larger than that of 
dislocations with non-zero components in both subsets, which we call dual-subset 
dislocations. Examples of the two types of dislocations are shown in Fig. [3] 
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Fig. 3: Two typical examples for the dislocations found using our procedure. 
In each counter-clockwise loop, the winding number rii is given by the number 
of red arrows, counting stripes crossed in the positive direction of b^*\ minus 
the number of green arrows, counting stripes in the negative direction, (a) A 
single-subset dislocation with Burgers vector (0,-1,0,-1) and (b) a dual-subset 
dislocation with Burgers vector (0,0,1,-1) (see text for definitions). 



To try and explain these observations, consider the free energy of the self- 
assembled crystal as an expansion in products of density mode amplitudes p(k) 
(see |Lifshit^|2QTT] ), 



^ip} = Yl E ^(ki,-..kn)p(ki 



•P(kn), 



(3) 



n=2ki...k„ 



where one can show that A{ki, . . . kn) = unless ki + . . . + kn = 0. We argue 
that products in the sum ([sj that contain high-intensity modes with large winding 
numbers have a greater contribution to increasing the free energy away from its 
minimum value in the perfect crystal. Accordingly, high-intensity modes tend to 
exhibit smaller winding numbers. Indeed, we find that all the winding numbers 
associated with the two brightest rings (see Fig. [ijb) and Fig. ^ are either or 
±1, whereas it is only on the 3^^ ring that we begin to see winding numbers that 
are either 0, ±1, or ±2. Moreover, owing to the linearity of the Burgers function, 
the fact that the winding numbers on the second brightest ring are at most of 
magnitude 1 prevents two adjacent peaks from the different subsets in the first 
ring from having non-zero winding numbers of the same sign. Because the ring of 
Bragg peaks, obtained by adding pairs of wave vectors separated by 60 degrees, is 
extremely weak [see Fig.[lJ^b)], there is no such constraint on the winding numbers 
belonging to the same subset of basis vectors. The fact that this constraint applies 
only to dual-subset dislocations, reduces their possible combinations and overall 
relative density. 
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(a) 



(b) 




Fig. 4: Inverse Fourier transform of a pair of Bragg peaks using filters of different 
size. For a circular filter with a radius of 25 pixels, used in (a), we find 12 disloca- 
tions, marked with red circles in (c). As we decrease the filter radius to 15 pixels, 
in (b), we find only 5 dislocations in the inverse Fourier transform in (d). 

A word of caution is in order regarding our approach for analyzing dislocations. 
Because the density of the dislocations is relatively high the Bragg peaks are not 
point-like but are rather spread as can be seen in Fig.jlfb). This means that some 
of the information about the dislocations may lie between Bragg peaks and may 
be lost if the filters are too small. Therefore, our approach is sensitive to the shape 
and size of the filter that we use around each Bragg peak. As we increase the 
filter size we obtain more information and potentially find more dislocations, as 
demonstrated in Fig. |4] We thus try to optimize the filter by gradually enlarging 
its size until the number of dislocations stops increasing substantially. 

Our approach for analyzing dislocations should be easily adapted to other 
systems even when the density of the dislocations is quite large, as one may expect 
for soft matter systems. Moreover, for dynamical systems that can be imaged in 
real time one can use our automated method to follow and quantitatively analyze 
the dynamics of the dislocations. 
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